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Abstract 

The phenomenology of the explicit CP violation in the Higgs sector of the next-to- minimal 
supersymmetric standard model (NMSSM) is investigated, with emphasis on the charged 
Higgs boson. The radiative corrections due to both quarks and scalar-quarks of the third 
generation are taken into account, and the negative result of the search for the Higgs bosons 
at CERN LEP2, with the discovery limit of 0.1 pb, is imposed as a constraint. It is found 
that there are parameter regions of the NMSSM where the lightest neutral Higgs boson may 
even be massless, without being detected at LEP2. This implies that the LEP2 data do not 
contradict the existence of a massless neutral Higgs boson in the NMSSM. For the charged 
Higgs boson, the radiative corrections to its mass may be negative in some parameter regions 
of the NMSSM. The phenomenological lower bound on the radiatively corrected mass of the 
charged Higgs boson is increased as the CP violation becomes maximal, i.e., as the CP 
violating phase becomes tt/2. At the maximal CP violation, its lower bound is about 110 
GeV for 5 < tan/3 ^ 40. The vacuum expectation value (VEV) of the neutral Higgs singlet 
is shown to be no smaller than 16 GeV for any parameter values of the NMSSM with explicit 
CP violation. This value of the lower limit is found to increase up to about 45 GeV as the 
ratio (tan/3) of the VEVs of the two Higgs doublets decreases to smaller values (~ 2). The 
discovery limit of the Higgs boson search at LEP2 is found to cover about a half of the 
kinematically allowed part of the whole parameter space of the NMSSM, and the portion is 
roughly stable against the CP violating phase. 
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I. INTRODUCTION 



The violation of CP symmetry has been with us for more than three decades since it was first 
observed in the neutral kaon system [1]. Within the context of the standard model (SM) [2] 
the CP violation is accommodated by virtue of a complex phase that exists in the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix for the charged weak current [3]. The origin of the complex 
phase has been remained in the standard model unexplained. A number of speculations have 
been proposed for possible explanations. One of the possibilities is that CP symmetry is violated 
in the Higgs sector, if there are at least two Higgs doublets [4]. It is known that super symmetric 
standard models require at least two Higgs doublets in their Higgs sectors. The two Higgs 
doublets give separately masses to the up-type quarks and the down-type quarks [5]. If the 
supersymmetry (SUSY) is exact in the model, where the SM particles and their superparticles 
have equal masses, CP violation would be absent in its Higgs sector, since certain restrictions 
from the SUSY are imposed on the Higgs sector. The supersymmetric model cannot however 
have an exact SUSY since the nature is not supersymmetric at electroweak level. If the SUSY 
is broken somehow in the model for phenomenological analyses, for example by introduction of 
the soft SUSY breaking terms, CP violation may arise in the Higgs sector of the supersymmetric 
standard model. 

In the minimal supersymmetric standard model (MSSM), many authors have investigated 
the possibilities of CP violation either explicitly or spontaneously. It is found that the tree- 
level Higgs potential of the MSSM is unable to break the CP symmetry in either way. Further 
investigations at 1-loop level of the MSSM show that the scenario of spontaneous CP violation 
is experimentally ruled out by the Higgs boson serch at the CERN e + e~ LEP, because it leads 
to a very light neutral Higgs boson [6]. However, at the 1-loop level, the scenario of explicit CP 
violation in the MSSM is found to be viable through the radiatively corrected Higgs potential 
[7]. Subsequently, a number of investigations have been devoted to the effects of explicit CP 
violation in the MSSM at 1-loop level to the charged Higgs boson [8] as well as to the neutral 
ones [9] by considering radiative corrections. It has been expected that at higher level the effect 
of explicit CP violation in the MSSM may lead to significant modifications to the tree level 
Higgs couplings to fermions and to gauge bosons. 

Meanwhile, the CP symmetry may also be broken if the Higgs sector is extended by intro- 
ducing at least one Higgs singlet. Then, the next-to-minimal supersymmetric standard model 
(NMSSM) is the right candidate for CP violation, which is basically different from the MSSM 
in that it has one extra Higgs singlet. For the spontaneous CP violation, it is well known that 
the tree-level Higgs potential of the NMSSM cannot provide the ground [10]. The situation is 
not improved at the 1-loop level by considering radiative corrections, if the top scalar quarks 
are degenerated in their masses [11]. In the NMSSM Higgs sector, spontaneous CP violation 
may only be realized by virtue of radiative corrections when the mass splitting effects between 
the top scalar-quarks are taken into account [12]. 

At the tree-level, unlike the MSSM, the scenario of explicit CP violation is acceptable by the 
NMSSM. By assuming the degeneracy of the top scalar-quark masses in the Higgs sector of the 
NMSSM, it is found that large explicit CP violation may be realized as the vacuum expectation 
value (VEV) of the neutral Higgs singlet approaches to the electroweak scale [13]. 

In this paper we investigate the effects of the explicit CP violation on the charged Higgs 
boson in the NMSSM. We calculate radiative corrections to the masses of both the charged 
Higgs boson and the neutral ones in the NMSSM due to the quarks and scalar-quarks of the 
third generation. The 1-loop effective potential which includes the effects of the mass splittings 
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in the scalar-quarks is employed. Assuming the nondegeneracy of the scalar-quark masses of 
the third generation in the 1-loop effective potential, a CP violating phase occurs from soft 
SUSY breaking terms, which penetrates into the radiatively corrected masses of the charged 
Higgs boson as well as of the neutral ones. Phenomenological implications and constraints by 
those CP violating Higgs masses are investigated over a wide region of the parameter space of 
the NMSSM, using the production processes of the neutral Higgs boson in e + e~ collisions. In 
particular, the results at LEP2 with y / i=200 GeV are used. 

Our paper is organized as follows. In the next section, we consider radiative corrections to 
the Higgs boson mass by using the 1-loop effective potential in explicit CP violation scenario. 
In the third section, we obtain phenomenological constraints on the charged Higgs boson mass, 
the VEV of the neutral Higgs singlet, by considering the neutral Higgs production at LEP2. In 
the last section, conclusions are discussed. 

II. HIGGS SECTOR 

The Higgs sector of the NMSSM contains two Higgs doublet superfields Hf = (H®, H^), Hj 
= (H2, H®), and a neutral Higgs singlet superfield N. Ignoring all quark and lepton Yukawa 
couplings except for those of the third generation quarks, the superpotential of the NMSSM 
including only terms with dimensionless couplings is written as [14] 

W = h t QH 2 t c R + h b QH x b c R + \NH X B 2 - \n* , 

where ht and hb are Yukawa coupling constants of the third generation quarks. The parameters 
A and k are dimensionless coupling constants. We define that H\H 2 = H^H^ — FL^ . The 
SU(2) doublet superfield Q T = (i^, bi) contains the left-handed quarks for the third generation. 
The SU(2) singlet superfields t c R and b c R are the charge conjugates of the right-handed quarks 
for the third generation. A global U(l) Peccei-Quinn symmetry would be explicitly broken by 
the presence of the cubic term of the Higgs singlet in the superpotential, which leads us to have 
a massless pseudoscalar Higgs boson at the tree level. 

At the tree level the Higgs potential of the model consist of F terms, D terms, and soft 
terms as 

V° = V F + V D + U soft , 

where 

V F = |A| 2 [(|#i| 2 + |# 2 | 2 )|Af + |#i# 2 | 2 ] + \k\ 2 \N\ 4 - (Xk* H1H2N* 2 + H.c.) , 

= l(3l + 9l)(\H 2 \ 2 - m 2 ) 2 , (1) 

VWt = m^J^ 1 | 2 + m^Ji7 2 | 2 + m 2 v |AT| 2 -(A^ Aj ff 1 F 2 Ar + i^7V 3 + H.c.) , 

with g\ and <? 2 being the U(l) and SU(2) gauge coupling constants, respectively. In the soft- 
terms, rrijj , m 2 H2 , and m 2 N are the soft SUSY breaking masses. A\ and are the trilinear 
soft SUSY breaking parameters with mass dimension. The role of the soft-terms in the Higgs 
potential is to break SUSY. The F terms come from the part of chiral auxiliary superfield while 
the D terms come from the part of gauge auxiliary superfield. 

In order to consider the explicit CP violation in the Higgs sector, A, A\, k, and A^ are 
assumed to be complex numbers. Among them, \A\ and kA^ can be adjusted to be real and 
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positive by redefining the overall phases of the superfields H1H2 and N. Then, we are left 
essentially with one phase: It may be chosen to be the one in Xk* = Xke 1 ^, which allows for 
explicit CP violation in the tree level Higgs potential of the NMSSM. We assume that the vacuum 
expectation values (VEVs) of three neutral Higgs fields H®, H®, and N are nonvanishing values, 
both real and positive, so that there is no spontaneous CP violation in the Higgs sector. 

The conditions that the Higgs potential becomes minimum at v\ = (H®), V2 = (H®), and 
x = (N) yield three constraints, which can eliminate the soft SUSY breaking masses mH 1 , mn 2 , 
and rriN from the potential V°. Note that although the phase <p is present in the potential, it 
disappears from the vacuum after applying the three minimum conditions for the VEVs. Thus, 
one does not impose a minimum condition on the Higgs potential with respect to the CP violating 
phase 4>. At the tree level, in consequence, the Higgs potential contains <f>, tan/3 = V2/V1, X, A\, 
k, Ah, and x as free parameters. 

The Higgs sector of the NMSSM has ten real degrees of freedom coming from two Higgs 
doublets and one Higgs singlet. After the electroweak symmetry breakdown takes place, three 
of the ten components correspond to a neutral Goldstone boson and a pair of charged Goldstone 
bosons, the other seven correspond to five neutral Higgs bosons and a pair of charged Higgs 
bosons. Transforming to a unitary gauge, one may express two Higgs doublets and one Higgs 
singlet as 

/ vi + Si + i sin (3 A \ 
\ sin/?C+* J ' 

/ cos/3C+ \ (2) 
V v 2 + S 2 + i cos (3 A ) ' 
{x + X + iY), 

where Si, S 2 , A, X, and Y are neutral fields and C + is the charged field. 

If CP symmetry be conserved in the Higgs potential, S±, S2, and X would become the 
CP even scalar Higgs fields whereas A and Y would become the CP odd pseudoscalar Higgs 
fields. The three orthogonal combinations among A and yield one neutral Goldstone boson 
and a pair of charged Goldstone bosons; they will eventually be absorbed into the longitudinal 
component of Z and W gauge bosons, respectively. We take the constraints = g 2 v 2 /2 and 
m z = (di + S'i)^ 2 /^ f° r the electroweak gauge bosons with v = \Jv\-\-v\ = 175 GeV. 

It has already been observed that the above tree level Higgs potential is significantly affected 
by radiative corrections due to the quarks and scalar-quarks of the third generation. Thus, 
we also take into account, in the presence of explicit CP violation, the radiative corrections. 
We employ here the method of the effective potential. The 1-loop effective potential including 
radiative corrections due to the quark and scalar-quark of the third generation is given by [15] 

where (i=l to 4) are the field dependent scalar-quark masses for the third generation, 

■Mq (q = b,t) are the field dependent quark masses for the third generation, and A is the 
renormalization scale in the modified minimal subtraction MS scheme. Note that the scalar- 
quark term is positive in the 1-loop effective potential whereas the quark term is negative. 
Because of the relative opposite sign between the scalar quarks and the quarks, if masses of the 
quarks and their corresponding scalar quarks are equal, then there would be no net contribution 
from the radiative corrections. However, this should not be the case since no light scalar particle 
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with mass comparable to light quarks has been observed: The quarks and scalar-quarks should 
not have equal masses. 

The top quark and bottom quark masses are given as m 2 = (htV2) 2 from hf^H^l 2 + |-£f^~| 2 ) 
and as m 2 = (h b v\) 2 from h 2 {\H^\ 2 + ji/^l 2 ) after the electroweak symmetry breakdown, re- 
spectively. By taking the basis {ti, t c R , &l, b R ) for the scalar-quark mass matrix of the third 
generation in Lagrangian density, the elements of the 4x4 Hermitian matrix, Mq, are explicitly 
given as 

M n = m 2 Q + h 2 \H Q 2 \ 2 + h 2 \H^\ 2 , 
M 22 = m 2 T + h 2 t \H%\ 2 + h 2 \H+\ 2 , 
M 33 = m 2 Q + h 2 \H«\ 2 + h 2 t \Hl\ 2 , 
M44 = m 2 B + h 2 b \H°\ 2 + h 2 b \H{\ 2 , 
Mn = h t (\NH° + A t H$*) , 

M 13 = -h 2 H°*H+-h 2 b H°Hr* , (3) 
M u = -h b \NH+ + h b A b H^* , 
M23 = h t XN*H^* - h t A t H+ , 
M 2i = h t h b H%H^* + h t h b Hf*H+ , 
M34 = -h b \NH$ - h b A b Hl* , 

where tuq, m 2 -,, and m 2 B are soft SUSY breaking masses and A q (q = t,b) are the trilinear soft 
SUSY breaking parameters. In the above expressions, the contribution from the D terms are 
not included. 

In order to accommodate the possibility of explicit CP violation at the 1-loop level, some of 
the coupling coefficients in the 1-loop effective potential are assumed to be complex. Note first 
that the complex phases of both ht and h b can later be absorbed away; they are thus assumed to 
be real. Next, we assume that tut = ms and A t = A b , on the ground that the leading radiative 
corrections come from the contribution of the top quark and top scalar-quark, although those of 
the bottom quark and bottom scalar quark for large tan (3 is not negligible. Accordingly, we may 
set the phase of At to be equal to that of A b . Consequently, at the 1-loop level, there remains 
only one relative phase between A t = A b and A. We name it as 4>i = 4>\ ~ 4>A q i and assume for 
simplicity that it is equal to the CP violating phase at tree level: 4>\ = 4>. 

The scalar-quark mass matrix for the third generation may conveniently be rewritten into 
four 2x2 block submatrices. Among them, the two off-diagonal block submatrices become zero 
matrices after the electroweak symmetry breakdown, since they do not contain the neutral Higgs 
fields. The two diagonal block submatrices may have non-zero mass eigenvalues, by substituting 
the VEVs for the neutral Higgs fields. Actually, by diagonalizing each of the block submatrices, 
one can obtain the top scalar-quark masses and m? 2 (m^ < m| 2 ) as 

m t + \( m2 Q + m T) T \l\( m Q ~ m T) 2 + m t ( A t + ^ 2x2 cot 2 (3 + 2A t \x cot (3 cos <p) , (4) 
and the bottom scalar-quark masses m| 3 and m| 4 (m| 3 < m~ 4 ) as 



m l + \ ( m Q + m V) T \j\( m Q ~ "4) 2 + m b ( A t + ^ 2x2 tan 2 /? + 2^Axtan/3cos</>) . (5) 

It is clearly the phase (p m the scalar-quark masses that is the possible source of the explicit 
CP violation. Notice that at the 1-loop level the phase (f> of the explicit CP violation would be 
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absent as long as the scalar-quarks of the third generation are degenerate in their masses. We 
assume that the left-handed and the right-handed scalar quark masses are not degenerate. 

The full Higgs potential of the NMSSM up to the 1-loop level may be decomposed into the 
neutral and the charged parts, according to which Higgs fields are involved. Then, the neutral 
part of the Higgs potential would yield the 5x5 mass matrix for the five neutral Higgs bosons, 
if it is differentiated twice with respect to the corresponding neutral Higgs fields. The charged 
part would also yield charged Higgs mass. 

Let us consider first the 5x5 mass matrix for the five neutral Higgs bosons. The elements 
of the matrix can be expressed, at the basis of (Si, 52, A, X, Y), as 

Mij = M% + bM\ 3 + bM% (i, j = 1 to 5) , (6) 

where come from the tree level Higgs potential, 5M^ from the top quark and top scalar- 
quark contributions, and 5M^ from the bottom quark and bottom scalar-quark contributions. 
We calculate them in an explicit manner to obtain the following results. First, the tree level 
Higgs potential gives 

M°! = (mz cos 0) 2 + Xx(A x + kx cos <p) tan (3 , 
M 22 = ( m z sin 0) 2 + Xx{A x + kx cos <fi) cot , 
M o = 2\x(A x + kxcos4>) 



sin 2/3 

M 4 ° 4 = (2kx) 2 - kxA k + ^v 2 A x sin 2/3 , 
Xv 2 

M 55 = T) — sin 2/3 + 3kxA k + 2Xkv 2 sin 2/3 cos <p , 

M° 2 = (A V - -m 2 z ) sin 2/3 - Xx(A x + kx cos <p) , 

M° 3 = 0, 

Mj° 4 = 2A 2 xv cos (3 — Xv sin 0(A X + 2kx cos (f) , 

M® 5 = — SXkvx sin f3 sin </> , 

M 2 ° 3 = 0, 

M 2 ° 4 = 2A 2 xv sin — Xv cos (3{A X + 2kx cos (f) , 

M25 = — 3>Xkvx cos (3 sin <p , 

M34 = A/cwxsin^ , 

M? 5 = Av(t4a — 2kx cos 0) , 

= 2A/cw 2 sin2/3sin0 , 

and the contributions of the top quark and top scalar-quark give 

3^A 2 x 2 A| ^m?, m|) 3 m 2 A^cos0 ? 
11 8tt 2 w 2 sin 2 [3 (m? - m? ) 2 87T 2 ?; 2 sin 2/3 7 1 *i ' *^ ' 

52 tl 

i 3m^4fA| g(m|, m|) 3m t 2 AxA f cos 2 2 

22 " STr 2 ^ 2 sin 2 /3 (m| - m? ) 2 + 167r 2 t> 2 sin 2 /3 tan (3 1 ^ ' ^ ' 



3mfA t A i2 bgjm|/m|) | 3m 4 ^ /mfm? 
47r 2 u 2 sin 2 /3 (m~ - m~ ) 87r 2 w 2 sin 2 (3 



j t '2' £11 j !H2 i oe 

T , 9 . 9 ,^ /_ 2 _ ™2 \ ^ n 9 9 • 9 A, iU & 
*2 51- 

,>iiil a j: , t, Me <^ 9{ m ^ 1 i n - t • aj .i./ mho 9 2 



3m 4 A 2 x 2 .4 2 sin 2 sC™?, ™|) 3m 2 Ax^cosv ,. - 
33 " SvrV sin 4 /3 (m| - m| ) 2 + 16vr 2 t; 2 sin 3 /? cos /3 ; ^ ' ) ' 
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mt _ 3m?A 2 A3 g(ml m|) 3rn|^cos0 2 , 
44 8vr 2 tan 2 (m| - m? ) 2 16vr 2 x tan /? 7 1 *i ' ^ ' 



3m 4 A 2 A 2 S in 2 ^gK, m|) gmg^cos^ , , 
8vr 2 tan 2 /3 (m? - m 2 ) 2 + 16vr 2 x tan /? 7 1 *i ' *' 

*2 tl 

3m 4 Ax^t A,- A,~ </(m? , m? ) 3m? Ax A cos . 



*2' 



t ■"»j^/i t a 1 -a 1 - 2 yy..^, ». fV ^ AXAf C QS g 2 2 , 

° IV1 V o T 2„2 c ; ln 2 / o / m 2 _ _2 \2 1 (W2„2 a : T1 2 „ ' "% > 



+ 



8tt 2 v 2 sin 2 f3 (m| - m| ) 2 16it 2 v 2 sin 2 /3 
3m 4 AxA fi log(m|/m|) 



87r 2 t; 2 sin 2 /3 (m~ — m~ ' 

A V *2 tl' 



3m^AVA f sin 0A fi g(m|, m| 2 ) 3m 2 Ax^sinc, , 2 

13 8vr 2 « 2 sin 3 (m? - m 2 ) 2 16it 2 v 2 sin /3 tan /3 1 K *V V 

t-2 tl' 

3m 4 A 2 xA| g(m|, m|) SmjA^ + Axcot/3) 2 

14 87r 2 «sin/3tan/3(m 2 - m 2 ) 2 16vr 2 i;sin/3 A«V "V> 

3m 4 A 2 xyl t sin0A t - i g(m|, m|) 3m 2 A^sin0 2 2 

15 8vr 2 t;sin/3tan/3 (m? - m? ) 2 ^ 16vr 2 t> sin/?tan/? 71 *V **> ' 
3m 4 AxA 2 sin^A t - 2 g(m|, m|) 3m 2 A x^sin^ , 2 

OiV/i 23 _ o_9..9 _:_3 o 7Z2 1^2 + 1C _2„.2„;„« /I 771 *,' m to> 



8TT 2 v 2 sm 3 3 (m~ — m 2 ) 2 IQit 2 v 2 sin/3 (l ' 
t 2 ti 



3m 4 AxAi sin <j> log(m|/m? , 
87r 2 w 2 sin 3 /3 (m~ — m? 



2 — m? " 
t 2 tr 



AMi Sm^A^A, g(m|, m|) 3 m 2 A^cos0 2 , 

24 8vr 2 v sin /? tan (3 (m| - m? ) 2 16vr 2 t; sin /? tan (3 7 1 *i ' J 
+ 3m 4 AA fi loglml/ml) 



87r 2 w sin /3 tan /3 (m~ — m? ) ' 

t 2 ii y 



3m 4 A^ 2 sin0A f - 2 gjm^ m|) 3m 2 A^sin0 2 2 

25 87r 2 «sin/3tan/? (m? - m? ) 2 167r 2 v sin tan /? A *i ' V 

t 2 ti' 

3m 4 AA t sin log(m| /m? ) 



87r 2 w sin (3 tan /3 



fm 2 — m? " 
v t 2 ti' 

,2 ™2 



, 3m 4 A 2 x^sin0A fi g(mj, rnjj 3m 2 AA t sin0 2 2 , 

34 " 8vr 2 t; sin 2 (3 tan (m? - m? ) 2 16tt 2 W sin 2 /3 7 lm *i ' ^ - 

3mf\ 2 xA 2 sin 2 5(^1 . ) 3m 2 AA t cos < 



8tt 2 v sin 2 /3 tan /? (m 2 - m 2 ) 2 16vr 2 t; sin 2 (3 fl ' * 2 ' 



3m^A 2 ^ t sin(/>Ar g(m| , 



^5 = - ' 2 ; o; "^2" -^9 , w 



87r 2 tan 2 /? ( — /;/ 



2 

ti ""ti 



and finally the contributions of the bottom quark and bottom scalar-quark give 

AM* 3m b A " A j g( m t , 3m 2 AxA f cos0 2 2 

11 8tt 2 « 2 cos 2 /? (m? 2 -m^) 2 167r 2 w 2 cos 2 /3cot /? 7 ^ ' "V 

+ 3m 4 A f A- b2 logCm^/m^) + 3m 4 / 

47r 2 v 2 cos 2 (3 (m~ — m~ ) 87r 2 u 2 cos 2 /3 °^ y m\ j 
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3m ' A x A ^ gK' m V 3mlXxA tC os<p 2 2 
8vr 2 t; 2 cos 2 /3 (to? - to? ) 2 8vr 2 ?; 2 sin 2/3 ^"V "V ' 

b 2 bl 

3to 4 A 2 x% 2 sin 2 4> to|) 3m 2 AxA f cos (f> 2 2 

8vr 2 ?; 2 cos 4 /3 (m? 2 - to^) 2 + 16ir 2 v 2 cos 3 /3 sin /? * ^ ' m &V ' 
3m 4 A 2 Ag i ff(m^ m|) 3 m 2 A^cos0 2 2 

8vr 2 cot 2 /3 (m? - to? ) 2 16vr 2 xcot/? A "V "V ' 
b 2 In' 



,2 ^2 

' 'dm£ A A t cos (p 2 2 

87r 2 cot 2 /3 ( m 2 2 - m 2 i ) 2+ 167r 2 xcot/? ^"V "V 



3m^ A 2 ^ 2 sin 2 <f> 9{m- > m s 2 ) 3m 2 A^ cos 

-/(mr m 



3m 4 AxA t A 5i A^ 5^, m? 2 ) 3m 2 Ax.4 t cos<K . 2 2 



87r 2 w 2 cos 2 /3 (to? — m? ) 2 167r 2 i> 2 cos 2 /3 61 ' 62 ' 
+ 3to 4 AxA^ \og{m\/ml) 



8ir 2 v 2 cos 2 /3 (to? — m? ) 

b 2 b 1 / 

3mi\xA 2 t sin <pA- b gimf m?J 3m 2 XxAt si n 2 

— — -U T ( 771- 111 ~ 1 

8vr 2 t; 2 cos 3 8 (m 2 -m?) 2 Wn 2 v 2 cos JK 61' 62 j 

b 2 bl 

_ 3m 4 Ax ^ sin ^(^/"() 
87T 2 v 2 cos 3 /3 (to? — to? ) ' 

02 bi 

SrojWV^ g^, m|) 3 m 2 A^cos0 2 ^ 2 
87r 2 vcos/3cot/3 (m? 2 -to^) 2 167r 2 t>cos/3cot/r l 
+ 3m 4 AA^ log^/fr^) 



87r 2 w cos /3 cot /3 (to? — m? x 

b 2 bi ' 



3m 4 A^ 2 sin^ 2 g(m^, to|) 3m 2 A^sin0 2 2 
87r 2 wcos/3cot/? (m? -to 2 ) 2 167r 2 vcos/3cot/3 Amb i' "V 

b 2 bi 

3m 4 A^sin0 ^(m^/m^) 
87r 2 v cos /3 cot 8 (to? — m~ ) 

b 2 bi 

3m 4 A 2 x 2 A t sin0Ar <7(m? m? ) 3m 2 Ax^sin0 „ 2 

- i i _ — _|_ _ t I tyi 

Sir 2 v 2 cos 3, 8 (to? — to?) 2 I6ir 2 v 2 cos /3 cot 8 61 ' 

b 2 bi 

3to 4 A 2 xA? 5 (m? to? ) 3to 2 A(A £i + Axtan^) 

-/(to, , 



87r 2 wcos/3cot /? (to? —TO?) 2 167T 2 UCOS/3 61 ' 62 

3m 4 A 2 x^ S in^ g(m? to? ) 3m 2 A^sin0 2 2 

— — -4- j 1 777,- 777- 

87r 2 u cos /3 cot 8 (to? -to?) 2 167r 2 vcos/3cot/3 b ^ b ? 

b 2 bi 

3m^A 2 xA t siii^A Si 9(.m 2 bi , m 2 ^) 3m 2 XA t sm<P 2 2 
8vr 2 «cos 2 /3cot/3 (to? - to? ) 2 16vr 2 t; cos 2 /3 ; ^1 ' "V ' 

b 2 bi 

3to 4 A 2 xA 2 sin 2 <f> djm^ , mfj 3m 2 XA t coscf) 2 2 
8vr 2 t;cos 2 /3cot/? (to? - to? ) 2 + 16tt 2 v cos 2 8 * ^ ' ' 

b 2 61 ' 



Sm^sin^ gjrr^, m|) 
^ " 87^ cot 2 /3 (m? - m? ) 2 " ( j 

In the above expressions, we have 

A j- = ^ cos + Ax cot /? , 

A t - 2 = A t + Ax cot (3 cos , (9) 
in <5M,*„- and 

A^ = At cos </> + Ax tan (3 , 

A- b2 = A t + Ax tan (3 cos 4> , (10) 
in bM\- , and we have the dimensionless functions / and g as 



f(mj, ml) 



m, — m 



2 2 
mi g A 2 " ~ m 2 ^2" 



+ 1 



mi + m? , m 2 



^m 2 ,™ 2 ) = ^l^ilog^f + 2. (11) 
mf — m\ 

Note that Mjg and are zero in M?- at the tree level. The former comes from the coupling 
between S± and A in the neutral Higgs potential, whereas the latter from the coupling between 
S2 and A. On the other hand, Aff 5 , M 2 ° 5 , M 3 ° 4 , and M 4 ° 5 are all non-zero. This implies that the 
scalar-pseudoscalar mixings at tree level may not occur between two Higgs doublets Hi and Hi- 
They can only occur either between the Higgs doublets and the singlet or in the cubic term of 
the Higgs singlet itself. At the 1-loop level, the scalar-pseudoscalar mixings between two Higgs 
doublets are generated by the radiative corrections due to the quark and scalar quark of the 
third generation. These scalar-pseudoscalar mixings are eventually responsible for the explicit 
CP violation in the radiatively corrected Higgs sector of the NMSSM. 

Every element of the mass matrix which describes the scalar-pseudoscalar mixings among 
the neutral Higgs boson masses are nonzero and proportional to sin <f> at the 1-loop level. Thus, 
the CP symmetry would be conserved in the Higgs sector if (f> is zero [16]. The maximal CP 
violation occurs when sin</> = 1. 

The five physical neutral Higgs bosons are defined as the mass eigenstates, obtained by 
diagonalizing the mass matrix at 1-loop level, by the help of an orthogonal transformation 
matrix. The elements of the orthogonal transformation matrix, Oij (i,j=l to 5), determine the 
couplings of the physical neutral Higgs bosons to the other states in the NMSSM. Let us denote 
the five physical neutral Higgs bosons as hi (i=l to 5). We assume that the mass of hi is smaller 
than hj if i < j. 

Note that when the vacuum expectation value (VEV) of the neutral Higgs singlet, x, is very 
large, one can derive the upper bound on the mass of the lightest neutral Higgs boson at tree 
level as [17] 

«,max) 2 = m 2 z + (AV - ml) sin 2 2/3 , (12) 

which is expressed in terms of only two parameters of the NMSSM: A and tan/3. None of A\, k, 
Ah, and x does enter in the expression. In particular, max does not contain the CP violating 
phase (p. 
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At 1-loop level the radiative corrections render the above expression. Assuming both explicit 
CP violation and large x, the upper bound on the radiatively corrected mass of the lightest 
neutral Higgs boson is given as 

(ml ? - K ?+ 3m ? (^cot/3A t - i+ AA- 2 ) 2 g 2 

+ 3m| (Axcot^+AtAf,) 1 /™| \ + _3mf_ Y r m l m L 
Air 2 v 2 (m~ — m 2 ) y m~ J 8ir 2 v 2 mf 



+ 



3mf (Axtan/?A 



8ir 2 v 2 (m~ 

02 




3mt (Axtan/?Ar + A t Ar ) /m^ \ 3m 4 m? m? 
+ 7-2V A * , 62 log + log(-^A ■ (13) 



47r 2 f 2 (m~ — m~ ) \ m~ / 8ir 2 v 2 mf 

v 62 61 \ 61/ b 

Note here that the upper bound is independent of the renormalization scale A. If <p = 0, there 
would be no CP violation, and taking cos 4> = 1 in the above expression for m\ x max would 
reduce to the one obtained with CP conservation. It is clear that the upper bound on the 
lightest neutral Higgs boson mass decreases as the CP violating phase (f> goes from to ir. 

We focus on the charged Higgs sector. In the NMSSM the CP violation cannot occur in an 
explicit manner in the charged Higgs sector because the pair of charged Higgs bosons are CP 
even. In a unitary gauge, the charged Higgs boson mass is obtained by differentiating the Higgs 
potential twice with respect to the charged Higgs field. At the tree level the charged Higgs boson 
mass in the NMSSM is given as 

( m ° c ) 2 = m 2 w - AV + X(A X + kx cos cp) <l + ® . (14) 

xan jj 

One can find that the relative sizes of the second and the third term would make m^ at the 
tree level either heavier or lighter than the W boson. Note that the dependence of m° c on the 
CP violating phase (/). If A\ be of the order of 1 TeV, m^ c would depend weakly on (j), and if A\ 
be as small as 0(0.1 TeV), it would strongly depend on (j). Moreover, unlike the neutral Higgs 
sector, the charged Higgs boson mass at the tree level depends on all the six parameters: tan/?, 
A, A\, k, x, and <f>. The dependence of m° c on A is not clearly visible in above expressions. 
However, it is easy to see that the mass of the charged Higgs boson would increase if any one of 
tan/3 (> 1), A\, k, x, and cos</> increases. 

The radiative corrections to the charged Higgs boson mass contain the contributions from 
the scalar-quark loops which depend on the scalar-quark masses of the third generation. Thus, 
one has to know the scalar-quark masses in advance. On the other hand, the scalar-quark masses 
depend on the charged Higgs field: The scalar-quark masses are given as the eigenvalues of a 
4x4 Hermitian mass matrix whose elements are coupled to the charged Higgs field. It is 
difficult to obtain an analytical expression for the scalar-quark masses in terms of the charged 
Higgs field. One can calculate the first derivative of the squared masses of the scalar quark with 
respect to the charged Higgs field. By replacing the result with VEVs, one can find a trivial 
fact. The only nontrivial part of the calculation is the evaluation of the second derivatives of 
the field dependent scalar-quark masses with respect to the charged Higgs field. 

The radiatively corrected mass of the charged Higgs boson is given as follows [18]: 

2 / \ 2 1 3m 2 m 2 2 2\ 

m c = (m v c ) + . /K, m b ) 

2tt z v z sin 2p 
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\ Bm\. + Cm\ + D 

v ,Wft "7 nK-<) ' (15) 

where m~. are the scalar-quark masses for the third generation. The expressions for B, C, and 
D are given as 

B = B^ + cos (3 sin j3{A t \x{ml sin 2 /3 + to 2 cos 2 /3) - m 2 m| sin 2/3} , 

C = C<t, + A 2 x 2 (m 4 sin 4 /3 + m 4 cos 4 /3) - XxA t m 2 (2ml + m 2 Q + mf.) cos 3 /3 sin/? 

- \xA t m 2 (2m 2 + mg + m 2 -) cos /3 sin 3 /3 

+ 2m 2 m 2 (m 2 + m 2 - 2,4 2 + mg + m| + A 2 x 2 ) cos 2 (3 sin 2 /? , 
I? = — A 2 x 2 m 4 (m 2 + m\) cos 4 /3 — A 2 x 2 m 4 (m 2 + m\) sin 4 /3 

- \xA t m 2 t (ml(A 2 + 2m 2 + A 2 x 2 ) - (m 2 + m 2 Q )(ml + m 2 -)) cos 3 /3sin/3 

- \xA t ml(m 2 t (Al + 2m 2 + A 2 x 2 ) - (m 2 + m 2 Q ){m 2 + to 2 -)) cos/? sin 3 /3 

- ra 2 TO 2 {2m 2 TO 2 + to 4 - - ^ 2 (m 2 + 4toq + to 2 + 2m 2 -) + (,4 2 + m 2 Q + A 2 x 2 ) 2 

+ (m 2 b + m 2 t ){m 2 Q + m 2 T + 2\ 2 x 2 )} cos 2 /3 sin 2 (3 , (16) 



where 



B^ = — A t Xx(m b sin /3 + m t cos (3) sin 2/3 sin — , 

= ^[4to 2 to 2 Ax + Ax(m 2 + to 2 )(toq + TOy) + Ax(m 2 — to 2 )(toq + my) cos 2/3 
+ AA t m 2 (3ml + to 2 + (to 2 - to 2 ) cos 2/3) (1 + cos^>) cos 3 /? sin/?] sin 2/? sin 2 ^ 



^(2Axto 2 (2to^t4 - 3A 2 m 2 - 2m 4 ) 

2 2\ 



+ 4to^Ax{toq(2to 2 + TOy) - TOj (to^ - 2to^ + 2X z x z )} 
+ Xx{A 2 mlm 2 + 4:(m 2 — m 2 )(3m 2 m 2 — TOqTO 2 -)} cos 2/3 
+ Atm 2 {32 \x Attn 2 cos 4 (3 sin 2 /3 cos 20 — 2\xAtm 2 cos 4/3 — XxAtm 2 cos 6/3 
+ (3to 4 — 3A 2 m 2 + 8to 2 to,q + 13to 2 to 2 + 3m 2 m^ + 5to 2 to^ — 5A 2 x 2 m 2 ) sin 2/3 
+ 2[m 2 (7m| + 3m 2 -) + to 2 (4toq + m 2 - - A 2 + to 2 - 3A 2 x 2 ) 
+ 4{m 2 (2m 2 + m 2 T ) + to 2 (toq + A 2 x 2 )} cos 2/3 
+ {m\{A 2 t - A 2 x 2 ) - (to 2 - to 2 )(to 2 + m^)} cos 4/3 
+ 16AxylfTO 2 cos 3 (3 sin /3] sin 2(3 cos 
+ 4(m 2 mg + 2m 2 m 2 + m 2 mj> + A 2 x 2 m 2 ) sin 4/3 

+ (^m 2 . - (mg - m 2 )(mg + m|) - AVm^sin 6/3}) sin 2/3 sin 2 | . (17) 

In the formula for the charged Higgs boson mass [Eq. (15)], the second and the third terms 
represent the radiative corrections: The second term comes from the quark loops of the third 
generation, and the third one from the scalar-quark loops of the third generation. As mentioned 
earlier, if the scalar-quark masses are degenerate in the 1-loop effective potential, radiative 
corrections to the charged Higgs boson mass would vanish and both the second and the third 
terms would be absent. 

III. PHENOMENOLOGICAL ANALYSIS 

We consider the neutral Higgs production in e + e~ collisions in order to investigate the effects 
of the explicit CP violation in the NMSSM, and to impose phenomenological constraints on 
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the parameter space using the LEP2 data. For the center of mass energy at LEP2, the main 
contributions to the production cross section for the neutral Higgs boson come from the Higgs- 
strahlung process and the Higgs-pair production process. The Higgs-strahlung process is 

(i) e+e~ -»• Z* -► Zhi (i = 1 to 5) , 

where the Z boson that is emitted together with the neutral Higgs boson in the final state is 
real for the center of mass energy at LEP2. When the next-to-lightest neutral Higgs boson mass 
is relatively small, the cross section of the Higgs-pair production process is comparable to that 
of the Higgs-strahlung one. The Higgs-pair production process is 

(ii) e + e~ — > Z* — > hihj for = 1 to 5, i ^ j) . 
The relevant couplings for the above two processes are 

G 2 ZZhi = ^Ij cos (30 il + sin f30 i2 ) 2 , 

G 2 Zhihj = — % r [{O j3 (cos pO i2 - sin 00a)- OaicospOjt- sin/30^)} 2 
J 4 cos uw 

+ {O i4 O j4 + O i5 O j5 } 2 ] , (18) 

where Ow is the Weinberg mixing angle, and Oij is the elements of the orthogonal transformation 
matrix that diagonalizes the mass matrix of the neutral Higgs bosons at 1-loop level. Thus, both 
couplings include radiative corrections and the effects of explicit CP violation. 

For numerical analysis, we input the values of relevant parameters in the NMSSM into our 
calculations. The masses of top quark [19] and bottom quark are fixed as 175 GeV and 4 
GeV, respectively. We assume that the mass of the lighter scalar-quark of the third generation 
is larger than the top quark mass. Then, we fix the renormalization scale in the effective 
potential at 500 GeV. We allow the remaining free parameters to vary within the following 
ranges: 2 ^ tan/3 ^ 40, < A ^ 0.87, and < k ^ 0.63. We assume that A\, A^, x, tuq, and 
wit may vary between and 1000 GeV, while A t between and 2000 GeV. The CP violating 
phase 4> is assumed to take any value between and tt. 

With those input values and ranges, we plot the contours of in Fig. 1 as functions of 
A x and A k , for tan/3 = 2, A = 0.12, k = 0.04, x = m Q = 1000 GeV, and m T = A t = 500 
GeV. We take the CP violation phase to be <j> = 7r/2, the maximal CP violation. In the same 
figure, the cross section at for scalar Higgs production in e + e~ collisions at y/s = 200 GeV is 
also shown, which is the sum of the cross section for Higgs-pair production process and that 
for Higgs-strahlung process. The hatched regions in the figure are where at is larger than 0.1 
pb and where my n ^ 0. If the values of the parameters of the NMSSM are such that they are 
represented by points in the hatched regions, it is likely that a scalar Higgs boson might be 
detected at LEP2 with ^fs = 200 GeV, if the discovery limit of LEP2 is assumed to be about 
0.1 pb. In other words, taking the results of the Higgs search at LEP2 into account, one should 
exclude the hatched regions in Fig. 1. 

On the other hand, the remaining region within the contour of m/ ll = is where at < 0.1 
pb. Therefore, if the parameters take the set of values that are represented by points in that 
region, a scalar Higgs boson might exist without being detected at LEP2 with ^=200 GeV. 
Consequently, the results of the Higgs search at LEP2 at -^5=200 GeV would not be able to 
put any constraints on mf ll . The LEP2 data does not exclude the possibility of the existence of 
a massless Higgs boson in the NMSSM with explicit CP violation. 
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In Fig. 2(a), we plot the charged Higgs mass as function of the CP violation phase in order 
to see the effects of the explicit CP violation in the NMSSM, for two different values of Ay. 200 
and 400 GeV. The remaining parameters are set as A^ = 50 GeV, tan/3 = 2, A = 0.12, k = 0.04, 
x = tuq = 1000 GeV, and tut = A t = 500 GeV. Those parameter values are consistent with the 
phenomenological constraints of LEP2 at \/s = 200 GeV. The two solid curves correspond to 
the tree-level masses of the charged Higgs boson, and the two dashed curves to the radiatively 
corrected ones. In both case, the upper one for A\ = 200 GeV and the lower one to A\ = 400 
GeV. We find that the tree-level mass of the charged Higgs boson is roughly stable against the 
variation of the CP violation phase while the 1-loop level one decreases slightly as as goes from 
zero to tt. From the figure, it is clearly shown that the radiative corrections to the charged Higgs 
boson mass are negative for ^ <fi ^ ^ an d that the amount of radiative corrections decrease as 
the CP violation phase goes from zero to tt. Another point to note is that for the whole range of 
<fi from to tt, the mass of the charged Higgs boson for A\ = 400 GeV is found to be larger than 
that for A\ = 200 GeV. This behavior does not change if the 1-loop contributions are included. 
Note that similar patterns are exhibited in the MSSM. 

In the lower part of the same figure, there are a dotted curve and a dash-dotted one. They 
are the plots of the radiatively corrected masses of the lightest neutral scalar Higgs boson, , 
as function of the CP violation phase, for the same set of parameter values. The dotted curve 
corresponds to for ^=200 GeV and the dot-dashed curve for ^4^=400 GeV. Unlike the 
charged Higgs mass, one can find that there occurs a crossover between the two curves as <fi 
moves from zero to tt: The mass for ^4a=200 GeV decreases as <fi increases, whereas the one for 
A\=400 GeV increases as <p increases. 

We repeat the same calcuations as Fig. 2(a) by changing the value of tan/3 from 2 to 10. 
The result is plotted in Fig. 2(b). In both figures, the values of the other parameters are same, 
and the figure captions are also the same. In Fig. 2(b), one can see that the behavior of the 
mass of the charged Higgs boson for tan/3 = 10 is quite similar to the one for tan/3 = 2, except 
that the mass of the charged Higgs boson for tan (5 = 10 is always larger by about 270 GeV than 
the one for tan [3 = 2, regardless of the value of <f>. From Figs. 2(a) and 2(b), we may infer that 
independently of the value of tan /3, the radiative corrections to the charged Higgs boson mass 
contribute more negatively as the CP violation phase increases. 

However, the behavior of for tan/3 = 10 is different from the one for tan/3 = 2. As we 
noted earlier, there is a crossover for tan/3 = 2 between m\ lx for A\ = 200 GeV and the one 
for A\ = 400 GeV as 4> moves from zero to tt. For tan (3 = 10, there occurs no crossover at all 
between the curves for A\ = 200 and 400 GeV. Figure 2(b) shows that m/ ll remains roughly 
constant as (f> varies, for tan/3 = 10. 

We perform the calculations repeatedly by changing the values of the relevant parameters 
to examine if peculiar behaviors might be found. For wide ranges of parameter values, no 
such peculiarity is found. The detailed numerical results are similiar to the behaviors shown in 
Figs. 2. In the presence of explicit CP violation, the radiative corrections to the mass of the 
charged Higgs boson is negative for ^ 4> ^ tt, whereas those to the mass of the lightest neutral 
Higgs boson may be either negative or positive contributions, depending on the particular set 
of parameter values. 

We turn to the Higgs production in e + e~ collisions. We feel that it is worthwhile to inves- 
tigate whether the LEP2 data at y/s = 200 GeV can put a phenomenological constraint on the 
values of the NMSSM parameters. In Fig. 3, we plot the lower bounds on x as function of tan/3 
for two different values of (j). The solid curve corresponds to (j) = (no CP violation) and the 
dashed one to <p = tt/2 (maximal CP violation). All the other parameters are allowed to vary 
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freely only if they satisfy the condition that at should be less than 0.1 pb, by using a random 
number generating function. 

The lower region of each curve is where the cross section for the Higgs production in e + e~ 
collisions at \fs = 200 GeV, a t , is larger than 0.1 pb, and thus should be excluded. We can 
say that x may by no means be smaller than 16 GeV for any set of values for other relevant 
parameters of the NMSSM, if the LEP2 data at yfs = 200 GeV is taken into account with a 
discovery limit of 0.1 pb. Thus, 16 GeV is a kind of universal lower bound on x in the NMSSM 
with explicit CP violation. Figure 3 also tells us that the lower bound on x increases up to as 
much as 37 GeV as tan (3 decreases to 10 down from 40. For 5 ^ tan j3 ^ 10, note that the lower 
bound on x for (f> = is larger than that for (p = it/ 2. 

We extend the preceding analysis to other relevant parameters of the NMSSM to see whether 
the LEP2 data at ^/s = 200 GeV can also put similar constraints on them. The procedure begins 
by generating random numbers for the parameter values, for a fixed value of the CP violating 
phase. For a partucluar set of parameter values, we calculate the masses of the neutral scalar 
Higgs bosons and check if the center of mass energy of the e + e~ collider, ^fs, is sufficient to 
produce the neutral Higgs boson via the aforementioned two processes. If the center of mass 
energy of the collider is larger than Et = mz + m\ li , then the Higgs-strahlung process is viable 
for the Higgs production. If ^/s is still larger than + mt lj , the Higgs-pair production process 
is also allowed. 

In this way we select the kinematically accpetable sets of the parameter values from 10 7 
points in the parameter space. The selected sets of the parameter values form a subspace within 
the whole parameter space of the NMSSM, which one may conveniently call as the kinematically 
accpetable parameter space of the NMSSM. The points in the kinematically acceptable space are 
further classified according to whether the cross section for the Higgs production they yield is 
smaller than 0.1 pb or not. If at is smaller than 0.1 pb, we finally choose the points. Otherwise, 
we discard them. The result is roughly an estimate of the fraction of the kinematically acceptable 
parameter space which can produce the neutral Higgs bosons without contradicting the LEP2 
data at ^fs = 200 GeV. We may call this subspace of the kinematically acceptable space as the 
LEP2-allowed parameter space of the NMSSM. 

In order to see the effect of the CP phase eft on the result, we repeat the procedure for 
different values of 4>, each time by generating 10 sets of the parameter values. The results are 
listed in Table 1. We see that about 44% of the kinematically allowed parameter space should be 
excluded. Thus, the size of the LEP2-allowed parameter space is approximately one half of the 
kinematically acceptable parameter space. Although the ratio of the excluded part in the case 
of the maximal CP violation (cp = ir/2) is larger than those of other values of 4>, the excluded 
fraction over the parameter space does not significantly depend on the CP phase cp. In other 
words, one may say that in the explicit CP violation scenario the data from LEP2 at ^fs = 200 
GeV can investigate about 44% for the kinematically acceptable NMSSM parameter space. 

With the knowledge on the allowed values of the parameters at hand, we proceed to estimate 
the mass of the charged Higgs boson in more detail. We generate randomly the values of the 
relevant parameters to select points in the LEP2-allowed space. In this way, 10 7 points are 
selected. For each point in the LEP2-allowed space, we calculate the mass of the charged Higgs 
boson, and plot it in the (mc, x)-plane. The results are shown in Fig. 4(a) for = (no CP 
violation), and in Fig. 4(b) for cp = ir/2 (maximal CP violation). Each figure contains the solid, 
dashed, dotted, and dot-dashed curves; they correspond to four different values of tan (3 = 2, 5, 
10, and 40. 
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Those curves are actually the boundaries of the plotted points in the (mc, x)-plane. Since 
there are no points plotted in the lower regions of those curves, the curves are clearly the lower 
bounds on the mass of the charged Higgs boson in the NMSSM. There are no curves for smaller 
values of x. This is because the smaller x is excluded, as can be seen in Fig. 3, by the condition 
of at < 0.1 pb, which is the assumed discovery limit for LEP2 at yfs = 200 GeV. 

In Fig. 4(a), we find that minimum values of the lower bound on the charged Higgs mass 
for tan/3 = 40, 10, and 5 are about 100, 95, and 85 GeV, respectively. If tan/? is as small as 2, 
the minimum of the lower bound on the charged Higgs mass might be almost zero. This implies 
that the NMSSM may even have a nearly massless charged Higgs boson, without contradicting 
the LEP2 data, if the relevant parameters have certain values. 

The minimum value of the lower bound on mc seems to increase sharply as tan/3 increases. 
In addition, the minimum value for tan/3 rises rapidly as x increases. For 300 ^ x ^ 1000 GeV, 
the charged Higgs boson might not have a mass smaller than about 110 GeV, regardless of the 
value of tan (3. The lower bounds for tan (3 = 5, 10 and 40 are smaller than that for tan (3 = 2 in 
the range of 300 s$ x s$ 1000 GeV. 

The above result is slightly different from, but essentially compatible with, the result of 
previous analyses [20]. In case of no CP violation, the previous analyses have estimated the 
lower bound on the charged Higgs boson mass is about 70 GeV. This number has been obtained 
under the assumption that the values of the relevant parameters are such that the mass difference 
between top scalar quarks are maximal. We do not assume to impose such constraints on the 
parameter values. The parameter values are randomly generated within the LEP2-allowed space. 
Thus, the top scalar quarks may have either maximally different masses or degenerate masses. 
In this sense, we may say that our result is more general. 

In Fig. 4(b), we repeat our analysis by substituting the value of 4> = by ir/2, to investigate 
the effect of maximal CP violation. Note that, as we assume that at the tree level the CP phase 
is equal to that at the 1-loop level, the CP violation arise maximally both at the tree level and 
at the 1-loop level. The behavior of the curves are basically the same as the case of <p = 0. we 
find that minimum values of the lower bound on the charged Higgs mass for 5 ^ tan [3 ^ 40 is 
about 110 GeV, and nearly stable if x is larger than about 100 GeV. The minimum values for 
(j) = tt/2 are larger than those for <p = 0, by about 25, 15, and 5 GeV for tan (3 = 5, 10, and 40, 
respectively. These increases may be regarded as the effect of the CP violation. 

If tan (3 is 2, and if x is about 75 GeV, there is a narrow possibility that the charged Higgs 
might be massless, without contradicting the LEP2 data, even if CP violation is maximal. Then, 
the minimum value of the lower bound on mc seems to increase abruptly as tan (3 increases from 
2 upward. For 300 ^ x ^ 1000 GeV, the charged Higgs boson might not have a mass smaller 
than about 110 GeV, regardless of the value of tan/3. 

IV. CONCLUSIONS 

In the next-to minimal supersymmetric standard model, we consider the explicit CP violation 
scenario in its Higgs sector and the subsequent phenomenology. We calculate the masses of both 
the charged Higgs boson and the neutral ones at 1-loop level, by taking into account the radiative 
corrections due to the quarks and scalar-quarks of the third generation within the explicit CP 
violation scenario. The effective potential at the 1-loop level, as well as the one at the tree 
level, is shown to have the CP violation phase. It would be absent if the top scalar-quarks are 
degenerate in mass. 
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In the whole parameter space of the NMSSM we establish a subspace with two constraints. 
The first constraint is that for the values of the relevant parameters in the subspace it is kinemat- 
ically possible to produce neutral Higgs bosons, and the second one is that the cross section for 
their productions does not exceed 0.1 pb. In this way, we establish the LEP2-allowed parameter 
space of the NMSSM, as we assume the discovery limit of LEP2 at y/s = 200 GeV. The lower 
bound on x, the VEV of the neutral Higgs singlet, is found to be about 16 GeV. It is shown to 
increase to 40 GeV as tan (5 approaches 2. 

We find that the LEP2 data do not exclude the possibility of a massless neutral Higgs boson. 
For the charged Higgs boson, we find that the radiative corrections to its mass are negative. For 
tan (3 = 2, we find that for certain values of x the lower bound on the radiatively corrected mass 
of the charged Higgs boson is as low as zero, whether CP is violated maximally or not at all. 
The minimum of the lower bound on the charged Higgs boson mass increases as tan/3 increases. 

For larger values of tan (3, the minimum value of the lower bound on the charged Higgs boson 
mass for maximal CP violation is found to be larger than the one for no CP violation. If x takes 
a value between 300 GeV and 1000 GeV, the lower bound on the mass of the charged Higgs 
boson in the NMSSM is about 110 GeV, for 2 ^ tan/3 ^ 40, whether CP is violated maximally 
or not at all. 
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TABLE CAPTION 



TABLE 1: Numerical results of at at \fs = 200 GeV for the parameter space of < tan f3 ^ 40, 
OA < 0.87, < k < 0.63, < A x , A k , x, m Q , m T ^ 1000 GeV, < A t < 200 GeV. The number 
of 10 6 points is kinematically allowed for the above parameter space at \fs = 200 GeV. 
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a t > 0.1 pb 
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553332 


553996 


545260 


567962 



FIGURE CAPTION 

FIG. 1. : The contours of and at in e + e~ collisions at y^s = 200 GeV, as functions A\ and 
A k , for 4> = tt/2, tan/3 = 2, A = 0.12, k = 0.04, x = m Q = 1000 GeV, and m T = A t = 500 GeV. 

FIG. 2. (a) : The parameters are set as A k = 50 GeV, tan/3 = 2, A = 0.12, k = 0.04, 
x = mg = 1000 GeV, and rar = At = 500 GeV. The dotted and dot-dashed curves are for , 
as a function 0, for v4>, = 200 and 400 GeV, respectively. The upper and lower solid (dashed) 
curves are for m° c (mc), as a function cf), for A\ = 200 and 400 GeV, respectively. 

FIG. 2. (b) : The same as Fig. 2(a), except for tan/3 = 10. 

FIG. 3 : The lower bound on x against tan/3 under the condition of a t < 0.1 pb for <p = 
and 4> = 7r /2- In this calculations, the remaining parameters are independently varied for 
< A < 0.87, < k < 0.63, < A x , A k ,m Q , m T 1000 GeV, and < A t < 2000 GeV. 

FIG. 4. (a) : The lower bound on mc under the condition of at < 0.1 pb as a function of x, for 
tan/3 = 2 (solid curve), 5 (dashed curve), 10 (dotted curve), and 40 (dot-dashed curve). In this 
calculation, the remaining parameters are independently varied for the same parameter space 
as Fig. 3. 

FIG. 4. (b) : The same as Fig. 4(a) except for <j) = tt/2. 
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FIG. 1.: The contours of and at in e + e collisions at yfs = 200 GeV, as functions A\ and 
A k , for 4> = tt/2, tan/3 = 2, A = 0.12, k = 0.04, x = m Q = 1000 GeV, and m T = A t = 500 GeV. 
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FIG. 2. (a): The parameters are set as A k = 50 GeV, tan/3 = 2, A = 0.12, k = 0.04, 
x = rriQ = 1000 GeV, and tut = A t = 500 GeV. The dotted and dot-dashed curves are for , 
as a function <j>, for A\ = 200 and 400 GeV, respectively. The upper and lower solid (dashed) 
curves are for m° c (mc), as a function <j>, for A\ = 200 and 400 GeV, respectively. 
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FIG. 2. (b): The same as Fig. 2(a), except for tan/3 = 10. 
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FIG. 3.: The lower bound on x against tan/3 under the condition of at < 0.1 pb for <j) = 
and (j> = tt/2. In this calculations, the remaining parameters are independently varied for 
< A < 0.87, < k < 0.63, < A x ,A k ,m Q ,m T < 1000 GeV, and < A t < 2000 GeV. 
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FIG. 4. (a): The lower bound on mc under the condition of at < 0.1 pb as a function of x, for 
tan/3 = 2 (solid curve), 5 (dashed curve), 10 (dotted curve), and 40 (dot-dashed curve). In this 
calculation, the remaining parameters are independently varied for the same parameter space 
as Fig. 3. 
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FIG. 4. (b): The same as Fig. 4(a) except for <p = tt/2. 
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